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We introduce a restricted four body problem in a 2+2 configuration ex- 
tending the classical circular Sitnikov problem to the circular double Sitnikov 
problem. Since the secondary bodies are moving on the same perpendicular 
line where evolve the primaries, almost every solution is a collision orbit. We 
extend the solutions beyond collisions with a symplectic regularization and 
study the set of energy surfaces that contain periodic orbits. 



1. Introduction 

One of the most important problems in celestial mechanics is the Sitnikov prob- 
lem [18], because this was the first restricted three body problem where the exis- 
tence of oscillatory movements was proved, as J. Chazy predicted in 1922 [3]. The 
Sitnikov problem is a generalization of the Macmillan problem introduced in [14] 
which is an integrablc problem and it has been studied by several mathematicians 
like Alcksecv [1] and Moser [15], Dankowicz and Holmes [5], Lacomba, Llibre and 
Perez-Chavela [12], Garcia and Perez-Chavela [7], among others. Some generaliza- 
tions of this problem include the Sitnikov problem in M 4 [12], the Sitnikov problem 
with three equal masses [6] , and recently the circular 4-body Sitnikov problem in a 
3+1 configuration [19]. In this project we study the 4-body Sitnikov problem in a 
2 + 2 configuration. In this configuration and for negative values of relative secon- 
daries' energy H < 0, in every solution the infinitesimal bodies collide. Therefore 
we consider collisions as elastic bouncing and we are interested in periodic solutions 
of this type, after applying the regularization process to continue solutions beyond 
collisions. 

Like other restricted problems, when the masses of infinitesimal bodies tend to 
zero, the system decouples and some singular terms vanish. This is the case we 
study in the present project. Instead of studying continuation of periodic orbits 
from circular to elliptic cases, we are interested in the conditions that must satisfy 
the values of fixed energy in order to accept resonant torus inside the hypersurface of 
constant energy. In a forthcoming work, we will study the transcendence conditions 
of the total fixed energy and its impact on the distribution of resonant tori. 

2. The 4-body Sitnikov problem 

The Sitnikov problem is a special case of the restricted three body problem where 
two massive bodies with masses mi = m-2 = h are evolving on keplerian orbits 
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around their center of masses, and there is an infinitesimal body that moves on the 
perpendicular straight line which passes across the center of masses of the massive 
bodies. The massive bodies are called primaries and the infinitesimal body is known 
as the secondary. The Sitnikov problem consists in determining the evolution of 
the body with infinitesimal mass under the attraction of primaries with Newtonian 
gravitational potential. 

The 4-body Sitnikov problem in 2+2 configuration (or double Sitnikov problem 
for short), is realized by the addition of one more secondary body on the perpen- 
dicular straight line where the first secondary evolves. In the general case, the 
secondaries have different masses = [i and = v with (i + 1 v and without loss 
of generality we can assume that v < fj, <C \- In this way, these bodies have no 
effect on the primaries' evolution, however with big positive masses of the secon- 
daries the dynamics of the system is very different to the circular classical Sitnikov 
problem. Of course, the two infinitesimal bodies interact between them under the 
Newtonian gravitational force. This is the subject of the dissertation work of the 
first author. 

The case with positive masses fj, > and v > will be called the reduced prob- 
lem, while the case with null masses fi = v = will be the restricted problem and 
it is called the 2+2 Sitnikov problem or alternatively the double Sitnikov problem. 
This work is related with the study of periodic orbits of the circular problem on 
the energy constant hypersurfaces. 

The potential of the reduced 2 + 2 problem in the general case is 

a v av 
V = , + — - H , 

Vgf + 1/4 yjq\ + 1/4 q 3 - 94 

the vector field is 



Mq 



dV 
dq 



and the Hamiltonian function is 



(1) K = -p'M-'p 



2 y/q* + 1/4 y/ql + 1/4 93-94 

(j, 



where p = (p3,Pi) are the conjugate momenta and M. = I ^ ^ j is the matrix 
of masses. 

In general, we can assign a correspondence rule to secondaries' masses in the 
form v = f(n) such that the lim^o f((A) = and then study the restricted problem 
which will depends on (i only. In our case we consider v = c [i with < c < 1. 

Remark 1. The case when 1 < c < oo is obtained by interchanging /i and v, 
so this analysis is valid for < c < oo. 

We obtain a new Hamiltonian function that now depends on c and /i as param- 
eters. 



1 Tu-l. a I 3 P 



(2) H = -p .M p -===== ===== — fi- 



2 v^s + 1/4 Vti + 1/4 93-94 

) 

/3 



where M. = ( ? „ ) , a = -A- and /3 = 1 
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FIGURE 1. The circular 4-body Sitnikov problem. 



At this point, we are interested in restating the problem from the symplectic 
point of view. We consider the open symplectic manifold (M, uj) defined as the 
cotangent bundle of the configuration space 1 Q = (M 2 \ A) where A = {q 3 = q^} is 
the set of singularities of H due to collisions. The manifold M carries the standard 
symplectic form 2 oj = ^ i dpi A dqi . We define the Hamiltonian system associated to 
the double circular Sitnikov problem, by H = (M,lu,Xh), where Xh is the vector 
field associated to the Hamiltonian function H : M — > R defined by (2). 

The Hamiltonian vector field Xh in local coordinates is as follows 

1 a q 3 (3 

93 = -P3, P3 = ; - Mt yo > 

■ 1 • £g4 , P 

94 = sPi, Pa = , + M7 



2.1. Regularization. To avoid the singularity in the Hamiltonian function and in 
the field Xh we extend analytically the equations to the hyperplanc q% = q±. We 
perform a symplectic regularization with the mapping p : M — > M defined through 
the generating function 

W(Q, p) = p 3 [qa + P^j + pa {Qa - cB± 
Then the mapping p : (Q, P) h-> (q, p) will be 

93 = Qa + P3 = nP 4 + -TT-, 

* H3 

(3) qA^Qi-a^, Pi=f 3P 4 -Il. 

* W3 

It is not difficult to show that p*(J2i^Pi A dqi) = J2i^Pi A dQi and, therefore 

p e Sp(M, lu) 



1 We say that this is the cophase space 

2 Some authors use lj' = £V dqiAdpi as the standard symplectic form. In the formal definition 
we consider the standar symplectic form as the exterior derivative of the canonical or Liouvillc 
1-form on M defined by 8 = YliPi^li- Consequently we have ui = d8 and ui 1 = —dS. 
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Also we consider the time rescaling 



(4) 



dt 
rf7 



aft Ql 



We will obtain a new function depending on the fixed value h = constant as a 
parameter in the following way: first we apply the change of coordinates defined 
by p(Q,P) = (q, p) and then the Hamiltonian function -ff(q, p) = H(p(Q,P)). 
Since p : M — > M is a symplectomorphism then the function in the new variables 
H(Q,P) = H(p(Q,P)) is again a Hamiltonian function. We fix the value of the 
function h = _ff(p(Q,P)) rearrange the terms and multiply by the rescaling time 
to obtain $L(H o p - h) = 0. 

The regularized Hamiltonian function is 

L = aft Qj (H - h) o p 

this Hamiltonian function depends on a, h as parameters and is valid only in the 
energy level L = for each h fixed. Specifically, if z = (Q3, Q4, P3, Pa) we apply 
the mapping p and the Hamiltonian function 



(H - h) o p(z) = 



2 {^ + af3Ql 



2a 



23 



(2Q 4 - aQlf 



(2Q 4 + (3Ql 



and after applying the time rescaling (4) we have 
1 



(5) 



{af3 PlQi 



Pi) 



2af3 z p 



-a/3Q 



2a 



2ft 



(2Q 4 



PQlf 



1 



(2Q 4 



aQlY 



1 



We write L/ l (z, fi) = L(z, n;h), and the dependence on the parameter a will be 
dropped because in this paper we only consider c = 1 as we will see below. 

We call to the triplet Ch = (M,Lo,X Lh ^ ^) the regularized system, where X^ h 
is the regularized Hamiltonian field. 

Although the form of the new Hamiltonian function is quite complicated, the 
advantage is that this function and the Hamiltonian vector field X^ h are regular on 
the boundary of M (specifically on the set A). Now we can obtain the limit when 
the mass p goes to zero 



lim L h {z,p) 
fjj — ^0 



and the effect is that the term ~2aft 2 p vanish. We can reverse the process, and 
since aftQ\ is not identically zero, we recover the Hamiltonian function in the 
original coordinates as follows 



(6) 



H 



1 . 

2 P M 



ft 



Vl'l + 1/4 + 1/4 
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Rewriting H_ = H(l + c) and considering the momenta pi = <ji for i = 3, 4, the 
original Hamiltonian function for the restricted case is 

(7) H = J_ \+c(U 



As we can see, the Hamiltonian function (7) corresponds to two uncoupled Sit- 
nikov problems. Figure 1 shows a diagram of it. 

The regularization permit us to continue analytically the solutions to the collision 
manifold q% = q±, however, if we want to study the problem as two uncoupled 
Sitnikov problems, we must give additional hypothesis in order to glue the solutions 
(contained in the same energy level) in a smooth way beyond the collision. When 
the secondaries have different positive masses < /1, < v = c/i with < c < 1, 
the elastic bouncing condition 

(8) 7)3 -V4, = -(v 3 - V4), 
and the conservation of linear momentum 

(9) P3 + Pi = P3 + Pi, 

implies the interchange of conjugate momenta at collision. (Here the terms pt and 
Vi = qi, i = 3,4 are respectively the momenta and velocities of the bodies after 
collision.) Therefore, for certain values of the masses fi and v there will be a 
discontinuity in the solutions of the original system (M, lu,Xh) for the restricted 
case. 

Using the equations (8) and (9) we can see the behavior of the system with 
Hamiltonian function (7) beyond collision. Writing expression (9) in the tangent 
space TqQ we have 

(10) CHV3 + f3v4 = ai>3 + f3v4, 

Let us solve for £4 in (8) and substitute in (10). Now, let us solve for £3 in the 
resulting equation to get 

(11) h — V3 - 2/3(v 3 - Vi), 
In the same way, we obtain 

(12) v 4 — v 4 + 2a(v 3 — V4). 

We write the system of two equations in matrix notation as v = Av where 
v = (v 3 , Vi) T and 



(13) A = 



1-2(3 2/3 
2a 1 -2a 



In the cotangent space T*Q this system is written as p = MAM~ 1 p. Making 
the computations we obtain AiAAi^ 1 = A T . Then the condition to continue the 
solutions in a smooth way beyond collisions using the transition conditions P3 = Pi, 

Pi = P3, h = Vi, and £> 4 = v 3 , is 

(i4) 

Thus, this is possible if and only if a = ft = 1/2 or equivalcntly [i~ v. We have 
proved the following 
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Figure 2. Solutions beyond collisions 



Proposition 2. In the circular double Sitnikov problem if // = v then the flow 
(pt(%) of the limiting case \x — > can be extended to a complete flow in a natural 
way considering crossing beyond collisions instead of elastic bouncing by means of 
the identification 

(15) q 3 - 94 >->• -(q 3 - q 4 ) 

which extends the Hamiltonian system to the whole phase space. 

On the other hand if we start with the conditions p± = p 3 and p 3 — Pi, and 
expanding p ~ A T p we arrive to 

Pi = (1 - 2/?)p 3 + 2ap 4 , 
p 3 = 2(3p 3 + (1 - 2a) P4 , 

from where 

(l-2 ( S)p 3 + (2a-l)p 4 = 0, 
and finally, since 1 — 2/3 = 2a — 1 we get 

(16) (2a-l)(p 3 +P4) = 0. 

This expression implies that if a ^ | (i.e. if /i ^ v) then the solution can be 
continued by changing signs p 3 = — p 3 and p 4 = —p^ only if p 3 + p 4 = 0. This 
situation is equivalent to reversing the solution after collision and this is the classical 
conception of elastic bouncing. In this case we will be interested in solutions that 
cross the 2-dimensional plane 

V = {q3 = 94} n {p 3 = -p 4 }, 

and this will be studied in a forthcoming paper. 

For convenience, we will analyze the problem in the original equations and this 
implies to use the expression (7) with c — 1. 

3. Action-angle coordinates and analytical solutions. 

Following Hoffer and Zehnder [10], we observe that if (M, o>) is an exact sym- 
plectic manifold of dimension 2n, there exists a 1-form such that ui = dX. For every 
symplcctomorphism ip <= Sp(M), the 1-form A — ip* A is closed and by the Poincarc 
lemma, locally there exists / : M ffi such that df = A — <p*\. 



ENERGY LEVELS OF PERIODIC SOLUTIONS OF THE CIRCULAR 2+2 SITNIKOV PROBLEM 



Integrating over a simple closed curve 7 we have J A = / <^*A, and finally we 



obtain 



A = 



A. 



[ A. An immediate 



We define the action on a simple closed curve 7 as ^(7) 

consequence of the above computations is that the action is invariant under sym 
plectomorphisms. 

This invariant property on simple closed curves permit us to construct a symplec- 
tomorphism for every periodic integrable Hamiltonian system (M,ui,H) that only 
depends on the values of the momentum map (although the term "action-angle" 
coordinates is generalized to non periodic Hamiltonian systems on non compact 
symplectic manifolds). 

Let H = h = /13 + /14 be the separable Hamiltonian function associated to the 
circular double Sitnikov problem. We can consider /13 < and /14 < and two 
closed orbits 7^ 3 , 7/i 4 associated to the relative energies and respectively. 
Then, there exist a symplectic change of coordinates <j) : M M where the trans- 
formed Hamiltonian function only depends on the action of the (simple closed) 
integral curves of each fundamental field. These coordinates are called action-angle 
coordinates, and are defined by 




Figure 3. Action-angle coordinates. 



2tt J lhi Q(hi) 
where fl(hi) = J^-. We have the following. 

Theorem 3. The action-angle coordinates for the double Sitnikov problem 
takes the form 

/2 

(17) J{hi) = -^(2E(ki)-K(ki)-IL(2klki)), 

(18) 0i{t;hi) = j£t(y,k)+0 o ,i, 
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where fli = 47r ^ 2 k' } ) (2E(ki) — K(ki) + II(2fc 2 , k t )) is the return time of the sec- 
ondaries, hi = ^ 2 + hi and 6o,i,i = 3,4, are constants determined by the initial 
conditions. 

Proof. The action is denned as the integral Ji(hi) = ^ § Pi dqi on a complete 
period. Since each hi is symmetric in pi and qt, we can integrate over a quarter of 
period and multiply it by four. 



(19) Ji(h) = 

I" JO 



hi H dq t 



1 



where q max is obtained when Pi = 0, then q max = y (-h t y 3- 

We construct a suitable change of variables (q,p) (y, s) that "normalizes" the 
integral, such that the following conditions hold: 

a) the Hamiltonian function takes the form h = y 2 /2 + a s 2 , 

b) for q = we require that s = 0, 

c) for q = q max we require that s = 1 and p = 0. 
The suitable change of variables 

(20) ___L== = (2 + -2, 



transforms the integrand of (19) in hi H — , = y/2 + hi\/l — s 2 . We write 

V V 9 i +1 / 4 

ki = y/2 + hi/2 and solve (20) for qi, then we compute dqi to obtain 



1 , „, 2 VI 



s 



2 



(21) , / H % = 9fcf , — ds. 



The integral (19) takes the form 



4V2, 2 I 1 VT^ds 



(22) Jif/ij) = -— k 

* Jo y/l - k 2 s 2 (1 - 2fc 2 S 2 ) 

This is a general complete elliptic integral. It is possible to write any elliptic integral 
in terms of algebraic rational functions of the independent variable, and the elliptic 
integrals of first, second and third kinds [9]. First we integrate by parts 



/ 



y/1 — s 2 ds s\/l — s 2 \/l — k 2 s 2 f s 2 Vl — fc 2 s 2 ds 



y/l - kfs 2 (1 - 2k 2 s 2 f l-2k 2 s 2 J y/l^{l-2k 2 s 2 )' 

We rewrite the last integral in the form 

1 f 2k 2 s 2 y/l - k 2 s 2 ds 1 „. , , If y/1 - k 2 s 2 ds 



2fc 2 J VT _ 7 2 (1 - 2fc 2 s 2 ) 2fc 2 v ' *' 2fc 2 J y/T^2(l - 2k 2 s 2 ) 

and again the last integral will be rewritten as 

1 f [1 + (1 - 2k 2 s 2 )]ds 1 1 2 

4F7 Vl-P,s 2 Vl- S 2 (l-2fc 2 S 2 ) _ 4^ F(Sl ' fcl)+ 4F n(2fc ' S *' fcl) - 
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Putting everything together we get 



VT^ds , Vl^sVl - k 2 s< 



y/l - k 2 S 2 (1 - 2fc2 S 2) 2 l-2fc 2 S 2 

(2E(s, kt) - ll(2kf, s , hi) - F(s, h)) . 

Finally, evaluating on the integration limits, we obtain (17). 

Now, the values of the period for each one of the secondaries is obtained in a 
straightforward way just by calculating the derivatives: 

(23) ^ = V 4(l-2fc2) ( 2E{h) K{h) + n(2 ^' ■ 

□ 

Note that the solution of the angle coordinates uses the time t = t(i>i,ki) that 
is not computed yet. This variable is obtained directly from the solution of the 
classical Sitnikov problem exposed by Belbruno, 0116 and Llibre in [2]. 

Theorem 4. The solutions for the circular double Sitnikov problem can be 
written as 

M ( h s(v 3 ) d(v 3 ) r- fc 4 s(v 4 ) d(v 4 ) r- 

where are functions of t obtained inverting the function 

( 24 ) * = 1 777-^3-7-^^ 



4(1 -2k 2 sn(v % ) 2 ) 2 

and s(vi) = sn(vi(t),ki), c{vi) = cn(fi(t),ki), d{vi) = dn(i>i{t) , hi) are the sine, 
cosine, and delta amplitude Jacobi elliptic functions, and ki = ^ 2 + hi for i = 3, 4. 

Proof. Since the solution of Newtonian Hamiltonian systems with one degree of 
freedom that have the form H = p 2 /2 — V(q) is 

r q dq 



(25) t + 1 

l qo ^2(h + V(q)) 
we can use the change of variables (20) in the differential 

, dq 1 ds 

(26) dt 



2^2 V / T~ r 7Vl - k 2 s 2 {\ ~ 2k 2 s 2 ) 2 ' 

yV+i/4 

and rescale the time for 4^dv to obtain 

dv 

, s dt 1 ds 
(27) — = — = and dv = 



dv 2V2{l-2k 2 s 2 ) 2 y/(l - s 2 )(l - k 2 s 2 )' 

We solve first the second integral directly because this is an elliptic integral of the 
first kind. The solution is s — sn{v, k) and this solution is substituted in the first 
integral to obtain the relation between the old and new times 

dv 



(28) 



2V2(l - 2k 2 sn{v, k) 2 ) 2 ' 
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Finally, solving for qi in (20) and substituting the solution for s = sn{v, k) we 
obtain 



s\/l - k 2 s 2 _ sn{v, k)dn{v, k) 
( ' Qi ~ 1 - 2k 2 s 2 1 - 2k 2 sn(u, k) 2 ' 

In order to obtain the conjugate momentum, we differentiate pi = ^ = ^7^- to 
get 

(30) Pi =^^ = 2V2kcn(v,k). 

dv at 

a 

It is possible to integrate the expression (24) with elliptic functions and elliptic 
integrals to obtain 



(3D* = 



1 - 2k 2 ) 



2E(v) - v + n(i/, 2k z ) - Ak 



2x Al2 sn{v)cn(v)dn{v) 



1 - 2k 2 sn(vf 



C. 



where C is an arbitrary constant of integration. In [4] the reader will find a nice 
and complete study of this function. 

From equation (31) the action-angle coordinates are completely described as 

J {hi) = —(2E(k l )-K(k l )-U(2k 2 ,k l )), 

7T 



7r / 2E{v) v + H>, 2k 2 ) - Ak 2 ^Z^lt^ 
2 I (2E(ki) - K(h) +U(2kf,ki)) 



7Q, i- 



From these expressions it is possible to deduce when Vi = K(ki) that 9i(K) = f + 
6> 0i i, therefore = AK(ki) implies 0i(4K) = 2tt + 8 i . Consequently, the solutions 
of the Hamiltonian subsystems generated by Hi = hi have period v = AK(ki) in 
the v variable and T(hi) = ^ in the t variable. This fact will be usefull in the 
analysis of periodic orbits for the circular double Sitnikov problem. 

Just one comment before to passing to the study of the level sets and periodic 
orbits of our problem. As the reader can observe, it is usually difficult to find 
the inverse transformation to have the new Hamiltonian function explicitly 
H( J, 6) = (7Jo</) _1 )( J, #), however since cf> : (q, p) >-> (J, 9) is a symplectomorphism, 

in particular dct (^-§f-^j = 1 an d applying the inverse function theorem, locally _1 
always exists. 



4. Hyper-surfaces of fixed energy 

In this section we describe the topology of the level sets for the Hamiltonian 
function (7) with c = 1, in terms of the momentum map and its image. Since 
for every point x £ Img(/i) its fiber ^i~ l {x) is a Lagrangian submanifold, we can 
decompose Img(/x) in smooth subsets with boundary to construct the foliations by 
hypersurfaces of constant energy. 
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4.1. Completely integrable Hamiltonian systems. A completely integrable 
Hamiltonian system is a Hamiltonian system (M,cj,Xh) and a set F = {Fi}f =1 
of first integrals for H = F± which are functionally independent and they are 
in involution (i.e., {Fi,Fj} = where {/,</} = lo{X f , X g ) denotes the Poisson 
bracket). In this case we call the set (M, u>, Xh, F) a completely integrable system 
in de sense of Liouvillc. 

In this context, the circular double Sitnikov problem is a completely integrable 
Hamiltonian system. Every first integral Fi, 2 < i < n generates an one-parameter 
family of symplectomorphisms by means of the exponential map. This one-parameter 
family can be realized as a Lie group G acting on the manifold M. 

We say that the action is symplectic if for every g G G, we have that the flow (p g is 
a symplectomorphism. Additionally, we say that the action is a Hamiltonian action 
if each of the fundamental fields is a Hamiltonian vector field. More specifically, 
the action tp t is Hamiltonian if there exists a map n : M — > Q*, from the symplectic 
manifold to the dual of the Lie algebra g = Lie(G) such that for every X £ g, the 
component \x x (p) := (fi(p),X) of p along X and for the fundamental vector field 
X* on M generated by the 1-parameter subgroup G° — {exp(tX)\t S R} C G, the 
relation 

dp x = i x »u 

holds, i.e., the function fi x is a Hamiltonian function for X* and fio ip g = Ad* o //, 
for all g £ G. 

Each fundamental field of an integrable Hamiltonian system is generated by one 
first integral Fi, such that {Fi, Fj} = for 1 < i < k. The application 

H = (H = F 1 ,...,F k ):M^g*=R k 

is called the momentum map and is defined from the symplectic manifold to the 
dual of the Lie algebra associated to the Lie group that acts on M. If k < n the 
system is partially integrable, however if k = n the system is Liouvillc integrable 
or completely integrable. 

We consider the momentum map p = (^3,^4) : (M,lu) — > M. 2 defined by 

(32) M = ^ pl ~7^W^ pl ~7^w)- 

In our case G has three possibilities 

1) G = S 1 x S 1 if ^(q, p) belongs to the third quadrant, 
11) G = S 1 x R if /x(q, p) belongs to the second or fourth quadrant, 
ni) G = R x R if /i(q, p) belongs to the first quadrant. 
In all three cases jj* = R 2 is obtained. 

It is well-known that the inverse image \i~ 1 (x) of each x G Img(/i) is a Lagrangian 
sub-manifold of {M,uS). If it is a compact set, it will be isomorphic to a torus. In 
other cases, it would be isomorphic to cylinders or planes according with the region 
where the point x lies (see Figure 4). 

In what follows we prove a result related to the image of hypersurfaces of con- 
stant energy of completely and separable integrable Hamiltonian systems under its 
momentum map. 

Lemma 5. Let (M, to) be an exact symplectic manifold of dimension 2n, and 
% = (M,ui, Xh) he a Hamiltonian system over (M,uj) with Hamiltonian vector 
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Figure 4. Values of the momentum map 



field defined by ix H u = dH. Suppose that there exists a symplectomorphism 
p : M — > M such that the new Hamiltonian function F = H o p is separable. 
Then there exists a Lagrangian Bbration tt : M — > R™ such that the hypersurfaces 
of constant energy H map to hyperplanes which are perpendicular to the vector 
1 = (1,1,..., 1) G R n . 

Proof. Since there exists p : M — > M, symplectomorphism such that F = Hop is 
separable then there exist global coordinates (Q,P) where F(Q,P) = H(p(Q,P)) 
separates in the form 

F(Q, P) = F 1 (Q 1 , Pi) + ■ • • + F n (Q n , P n ), 

and Fi(Qi,Pi) = constant are n first integrals for Xp. Moreover {F,Fi} = 
for i = l,...,n. The Hamiltonian system is integrable by quadratures and we 
can consider the combined flow cp* of all the Hamiltonian vector fields Xp i as a 
Hamiltonian action of the Lie group G = M. k x f n ~ k on (M, lj) for some 1 < k < n. 
The Hamiltonian action G x M — > M induces a momentum map 

H = {F u F 2 ,...,F n ):M->g*, 

where g* = R" is the dual of the Lie algebra associated to G = tp*. Its image 
Img(/i) C R" is a convex polyhedron or cone whose vertices are the extremal values 
of p as was studied by Guillcmin and Stcnbcrg in [8]. The image /x(S) of every 
regular hypersurface of constant energy S/j = under the momentum map 

is a convex subset ICh of the linear affinc subspacc of codimension 1 of q* 

x\ + X2 H 1- x n — h = 0, 

where x = (xi, . . . , x n ) € R™. We can write (x, 1) — h = 0, and in particular we 
have 

JC h := {x G R" |x G (Img(M) PI {(x, 1) = h})} C {x G R"|(x, 1) = h}. 
Then ir := p o p^ 1 : M — > R™ is the smooth map we are looking for. 
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ii 



g* 

Finally, we know that the fibers £ = p^ 1 (x) with are Lagrangian subman- 

ifolds of M for every x € lmg(p), this implies that uj\s, = 0. Since p^ 1 € Sp(M,ui) 
then 

(33) u(x,y)=u{p-\x),p- 1 (y)) = 0, V x,y € £ 

therefore £ = p _1 (£) is a Lagrangian submanifold. We conclude that ir : M — > M™ 
is also a Lagrangian fibration. □ 

Remark 6. It is important to note that the interior points of the set fCh 
correspond to Lagrangian submanifolds of M. On the other hand, the points lying 
on the boundary dlCh correspond to isotropic submanifolds that we can think as 
degenerate Lagrangian submanifolds. We mean that if a: € dJC then Lj\^-ir x ) = 0. 
The isotropic submanifolds have the form R r x T s with < r, s < n and r + s < n 

Corollary 7. The hypcrsurfaces of constant energy of the circular double 
Sitnikov problem under the momentum map (32) correspond to segments of lines 
with slope m = — 1 in lmg(p) C R 2 . 

Proof. Since the circular double Sitnikov problem is a Hamiltonian system de- 
fined on M = R 4 , the image of the momentum map (32) is a subset of R 2 and 
the affine subspaces perpendicular to 1 R 2 = (1,1) are the straight lines with slope 
m = — 1. □ 

4.2. Level sets of fixed energy. In order to describe the surfaces of constant 
energy and their foliations, we consider the separable Hamiltonian function in the 
following form 

H = h = /13 + hi, 

where hi corresponds to an energy level of the circular classical Sitnikov problem, 
for each i = 3, 4. Then we consider the image of the momentum map (32) and finally 
we construct the foliation following the straight line associated to each surface of 
constant energy in Img(/i). 

From the solutions for the classical Sitnikov problem [2], we know that hi, for i = 
3, 4, is defined in [—2, 00) and the orbits have the following behavior: if —2 < hi < 
the circular Sitnikov problem has periodic orbits, for hi = it has a parabolic orbit 
and for hi > it has hyperbolic orbits. Due to the restriction on each relative 
energy hi, for i = 3,4, the total energy h = hs + hi has the image (—4, 00) and we 
have the following topology (see Figure 5): 

• If h = —4 this level does not exist in the real problem q\ = because the 
secondaries are in the same place at the same time (impossible). 

• If —4 < h < —2 the energy levels corresponds to spheres S 3 foliated by tori 
T 2 and two singular closed curves. 

• If h = —2 the energy surface is a 3-sphere without four points. 
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FIGURE 5. Lines associated to surfaces of constant energy 

• If— 2</i<0 the energy surfaces are 3-spheres with four discs D 2 as 
boundaries. 

• If h = the foliation contains two disjoint cylinders with four planes in the 
middle point (when h% = h± = 0). 

• If h > the foliations contains cylinders and planes. 

A few of these foliations are shown in figure 5. 

The most interesting energy levels are when h = — 2 and h = because these 
are bifurcation values for the topology of the constant energy surfaces. Other 
interesting energy levels are —4 < h < —2 because the energy surfaces are 3- 
spheres foliated by 2-tori and they have all the solutions bounded, and this gives 
the possibility of finding interesting periodic solutions that will be preserved under 
small perturbations of the eccentricity e for the kcplerian solutions of the primary 
orbits, or perturbations on the mass parameter /i of the secondaries. 

5. Periodic solutions for the circular double Sitnikov problem 

At this point we have shown that every solution of the 2+2 Sitnikov problem 
has the form 



with v = v{t). When the values of the momentum map are in the third quadrant 
the evolution of the system is bounded. In this region it is possible to have peri- 
odic orbits of the four bodies under specific conditions. In what follows we give 
some definitions and we establish the conditions that produce periodic orbits in the 
circular double Sitnikov problem. 

Definition 8. We say that ip(t) is a periodic solution of period r with r > 
if ip(t + t) = <p(t) for all t S M and there does not exist f <G (0, r) such that 
ip(t + f ) = f(t), i.e., r is the minimum period. 

Since the solutions of the double Sitnikov problem are in terms of the Jacobian 
elliptic functions sn{v, k), cn(v, k) and dn(v, k) which are defined on the Riemann 



(34) 
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surface generated by two primitive periods 2K and 2iK' in general, they accept 
complex arguments and modules. In fact, these functions are analytic functions 
in the module k € (C \ { — 1,1}), but just if k £ (R \ {1,1}) its image is real. 
Here, K = K(k) is the complete elliptic integral of first type. Therefore, the 
body with position qi will have a return time v = AK(ki) in the rescaled time and 
t(4:K) = T(hj) in the real time t = t(v), for i = 3, 4, 

(35) T(h) = 2(1 _^ 2fc2) [2£(fc)-A-(fc)+II(2fc 2 ,fc)], 

where k = ^ 2 + h and K, E, and II are the complete elliptic integrals of first, second 
and third type respectively. 

We will need some more properties about the function T(k) = T(h), which are 
summarized in the following result. 

Theorem 9 ([2]). Let T(k) = T(h) be the period of the solution of the circular 
Sitnikov problem with energy h; then the following statements hold. 

(1) Hm^_ 2+ T(/ l ) = ^. 

(2) lini/j^o- T{h) = oo. 

(3) §>0, W*e(-2,0). 

(4) hm^„ 2+ § = 

(5) lim^o- ^ = oo. 

The proof of this theorem follows directly from the definition of the period 
T(k(h)) = T(h) as function of h. We refer the reader to [2] for details. 

With these elements, we will characterize the periodic orbits of the double Sit- 
nikov problem. We will use the notation (p, q, n) = 1 to mean that the greatest 
common divisor is gcd(p, q,n) = 1, in other words, that the three numbers have not 
common factors at the same time. 

Proposition 10. For every periodic solution of the double Sitnikov problem 
there exist 3-plcts (p, q, n) G Z 3 such that (p, q,n) — 1, and p > and p > 
holds. The periods of these solutions are related to the partial energies by 

r = 2pir = qT{hs) = nT(h4). 

Remark 11. The couples (p, q) and (p, n) are not necessarily coprime, however, 
at least one of the three combinations (p, q), (p, n), (q, n) must be coprime to assure 
that (p, q,n) = 1. 

Remark 12. The energy surface of the double Sitnikov problem that accepts 
periodic solutions with period r = 2tt = T{h^,) = T(hi) is a non compact hyper- 
surface. The value of T(h) in h = —1 is 

- V2^(i)-F(i) + n(ii) 

The numerical estimation of this value is 

T (~ 1 ) = 0.824429907123718 < 1. 
2tt 
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Since the function T(h) is an increasing function of h thus T(hi) = 2ir is obtained 
for hi > —1 and — 2 < h.3 + hi, therefore E = i? _1 (/i3 + hi) is not compact (see 
Figure 5). 

Definition 13. We say that an energy surface E = i? _1 (h) accepts a periodic 
solution if there exists p,q,n € N with the following properties: 

PI. (p,q,n) = 1, 
such that 

E = H- 1 (V 1 fevr) +T- 1 (^2tt 

We will write E/j = iJ _1 (h) in order to make clear the dependence on h. 
We denote the set of fixed energy surfaces that accept periodic orbits as 

Tl = |e = H-\h*)\K = T- 1 [j^ 1 ^ (^ 27r ) ^^-P 2 holds 

Theorem 14 ([11]). In the circular double Sitnikov problem there exists a 
countable number of energy surfaces E £ 9H that contains resonant tori foliated by 
periodic orbits. Moreover, the set of values h* 6 H(M) C K such that E^, € 2Tt is 
dense in (-4,0) and have zero measure in K. 

It is a well-known result that resonant tori form a dense set in the image of the 
momentum map. However, Pugh and Robinson proved in 1983 [16] that generically 
the periodic orbits of Hamiltonian systems are dense in any open set contained in 
the union of compact and regular energy surfaces. Moreover, they argued that 
using a Fubini's argument, this result apply for any given compact and regular 
energy surface. In contrast, last theorem assures that there exists a set of values 
h* € M of full measure such that E/ u ^ 9Jt. That is a generic behavior of completely 
integrable Hamiltonian systems. 

The proof of Theorem 14 is an immediate consequence of the following two 
lemmas that we now state and prove. 

Lemma 15. For each n £ N the circular double Sitnikov problem has periodic 
solutions of period 2nir. 

We will just exhibit at least one periodic solution of period t = 2Nir. This is 
immediate from the fact that there exists such periodic solutions in the circular 
(classical) Sitnikov problem. 

Proof. For any iV £ N we can choose the combination p = N and q = n = 1 
that produce 

(p,q) = l and (p,n) = l, 

with 
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and Proposition 2.8 in [4] assures that there exists h±,h2 £ (—2, 0) such that 

T(h 3 )= 2 -^ and T(h 4 ) = ^. 

q n 

Then the hypersurface H (ha+h^) contains a torus foliated by a family of periodic 
orbits with period 

t = 2nN = T{h 3 ) = T{hi). 

□ 

The following lemma is about the finiteness of resonant tori foliated by periodic 
orbits of prescribed period t. 

Definition 16. We define the totient function or Euler's phi function <p(p) of 
an integer p by 

P(P) = ^Ilf 1 ^) 

n\p V 7 

where the product runs on all n coprime to p. It represents the number of positive 
integers less than or equal to p that are coprime to p. 

Lemma 17. For each N £ N fixed, the circular double Sitnikov problem have a 
finite number of tori foliated by periodic orbits with period r = 2Nir. The number 

8Ncp(N) + 

q < 2y/2N, 

is an upper bound (although is not an optimal bound). 

Proof. For each p £ N fixed there exist 3-plets (p,q,n) £ N 3 , where properties 
PI and P2 of Definition 2 holds. Therefore, we search for the number C p of 3-plets 
(j>,q,n) = 1 coprimes. It is easy to see that for every q < 2y/2p and (p, n) = 
1, the 3-plet (p,q,n) does not have common divisors. These triplets are exactly 

(2V2 P ) ■ (2V2^(p)) = 8Mp). 

Additionally, we must add all the couples (q, n) coprime such that (p, q) and 
(p, n) are not coprime. This means that for each integer q < 2\f2p with (p, q) ^ 1 
we must add the number of coprimes <p(q). Then we have 

(37) C p <8p<p(p)+ 

q < 2^2p 
(P,3)^l 

Finally we must eliminate the elements that are in both sets, however the number 
(37) is an upper bound of the 3-plets (p,q,n) £ N 3 where properties PI and P2 
hold. 

The 3-plet (p,q,ri) £ N 3 induces a point x = (2ir^,2ir^) £ (T(/i 3 ), TQia)) such 
that the Lagrangian torus T = oT^ 1 ){x) is foliated by periodic orbits of period 
2Nir, therefore it is a resonant torus Tn es C M. □ 
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Proof of Theorem 14- The first part of the theorem is a consequence of the fact 
that the countable union of finite sets is a countable set. Using Lemmas 2 and 3 we 
have that the number of resonant tori are countable, and since each torus belongs 
to exactly one energy surface, the set 9Jt is countable too. 

Now we must to prove that the set of values /i* of energy surfaces with resonant 
torus is dense in (—4, 0), and have zero measure there. We define the map T : Q* — > 



For each rational point y G Img(T) with y = (~, ~), (r, s) = 1 and (u,v) = 1, 
we construct the point (^, ^t, y ) <G N 3 where g = gcd(ru, su, rv). Since this point 
fulfills properties PI and P2 of definition 2, there exists a resonant torus foliated 
by periodic orbits with period 



The set of rational values of T defined by RP := Img(T) H <Q> 2 is a dense subset of 
zero measure in Img(T). The mapping T is continuous and then T^ 1 {RP) C g* 
is a dense subset in the image of the momentum map p. Now we construct the 
function H : g* — >• M such that sends x = (ft.3,/14) 1 — >■ + /14. It is immediate 
that "H(T -1 (RP)) C (—4, 0) is a dense subset by continuity, and have zero measure 
since RP is a countable set. □ 



In this section we use some facts about the transcendental number theory related 
to the transcendence of the periods of elliptic functions, in order to characterize the 
values li, el such that we have € 9Jt. 

The results of the last section can be restated as follows 

Theorem 18. Every point x <S Img(/i) is the projection of a resonant torus 
foliated by periodic orbits of the circular double Sitnikov problem if, and only if 
T(x) is a rational point. 

Proof. Suppose that % is a resonant torus, and /x(T) = x G q* with x = (/13, hi). 
It means that there exists a number r € M + such that <j(t + r) = <j(t) holds for 
every solutions on 1. Moreover, r = 2nir for some n S N. Since t is the (minimum) 
period, then there exists p, q € N such that pT(hs) = t and qT(hi) = r with 
(n,p,q) — 1. We obtain that T(h 3 ,h 4 ) — (n/p,n/q) £ Q 2 is a rational point. □ 

Now, we want characterize the values of the relative energies € (—2,0) 

which produce resonant tori. Using some relations between the elliptic integrals 
and functions of Jacobi we have the following expression for the complete elliptic 
integral of third kind 

n(A'(fc),2fc 2 ,fc) = K(k)E(2k 2 ,k)~2k 2 E(k) 

(formulae (3.8.32) and (3.6.1) in [13]). Therefore, from (35) the condition for T(x) 
be a rational point is equivalent to 



by 





6. A CONJECTURE 



(38) 



Ti 1 



[2(K) 2 E(k l ) - (1 - E(2k 2 , fc i ))A'(fe i )] e Q, 



2tt 4ttv / 2(-/i 1 ) 
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where E(ki),i = 3,4 is the complete elliptic function of second kind, (k^) 2 = 1 — kf 
is the complementary modulus, and E(2kf,ki) is the incomplete elliptic function 
of second type with argument 2kf and modulus fcj. In the last formula, the ratio 
Ti/(2ir) is expressed in terms of elliptic functions of first and second kind only. 
Then, it is possible to apply some results on transcendental number theory due to 
Schneider [17] in order to characterize the values of hi and ki such that expresion 
(38) holds. 

Conjecture 19. If the circular double Sitnikov problem has a periodic solution 
with period r = 2nir, n £ N then the relative energy values belongs to the field Q(k*) 
where k* G {R\ A}. This field is an extension with degree of transcendence 1 over 

Q. 

It means that all the constant energy values where the resonant tori lie are 
algebraically dependent. 
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